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ABSTRACT 

We investigate the sensitivity of epidemic behavior to a 
bounded susceptibility constraint - susceptible nodes are in¬ 
fected by their neighbors via the regular SI/SIS dynamics, 
but subject to a cap on the infection rate. Such a constraint 
is motivated by modern social networks, wherein messages 
are broadcast to all neighbors, but attention spans are lim¬ 
ited. Bounded susceptibility also arises in distributed com¬ 
puting applications with download bandwidth constraints, 
and in human epidemics under quarantine policies. 

Network epidemics have been extensively studied in lit¬ 
erature; prior work characterizes the graph structures re¬ 
quired to ensure fast spreading under the SI dynamics, and 
long lifetime under the SIS dynamics. In particular, these 
conditions turn out to be meaningful for two classes of net¬ 
works of practical relevance - dense, uniform (i.e., clique¬ 
like) graphs, and sparse, structured (i.e., star-like) graphs. 
We show that bounded susceptibility has a surprising im¬ 
pact on epidemic behavior in these graph families. For the 
SI dynamics, bounded susceptibility has no effect on star- 
like networks, but dramatically alters the spreading time 
in clique-like networks. In contrast, for the SIS dynamics, 
clique-like networks are unaffected, but star-like networks 
exhibit a sharp change in extinction times under bounded 
susceptibility. 

Our hndings are useful for the design of disease-resistant 
networks and infrastructure networks. More generally, they 
show that results for existing epidemic models are sensitive 
to modeling assumptions in non-intuitive ways, and suggest 
caution in directly using these as guidelines for real systems. 
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1. INTRODUCTION 

Epidemic processes - stochastic models for spreading in 
graphs - play an important role in a variety of domains 
ranging from social networks (opinion spread), epidemiol¬ 
ogy (disease spread) and distributed computing (informa¬ 
tion spread). Though there are several models in literature 
for such processes, a characteristic feature of these is the 
spread of the process from ‘infected’ to ‘non-infected’ nodes 
along the edges of the network. In this work, we investi¬ 
gate the behavior of such processes under the assumption 
that the maximum susceptibility of a node to infection from 
its neighbors is bounded. This constraint seems natural (and 


we present several settings where it is so) - however we show 
that imposing it brings to light some surprising behavior of 
epidemic processes in different networks. 

In this paper, we focus on two major classes of epidemic 
processes: the SI process where nodes never recover from 
infection, and the SIS process where nodes recover, but can 
get reinfected. The SI or ‘Susceptible-Infected’ dynamics [1, 
2, 3, 4] is the most common model for one-way dissemina¬ 
tion. Nodes exist in 2 states - ‘susceptible’ (S) and ‘infected’ 
(/). Infected nodes never recover, and try to infect each 
neighboring susceptible node at a rate P (after a random 
delay, drawn from an i.i.d. Exponential{P) distribution). 
We refer to p as the virulence of the infection process. The 
SI dynamics terminate when all nodes are infected, and the 
corresponding time interval is called the spreading time. In 
particular, we are interested in networks with small spread¬ 
ing time (ideally 0(1) in the number of nodes). 

A small modification of this model gives the popular SIS 
or ‘Susceptible-Infected-Susceptible’ dynamics [5, 6, 7, 8, 9, 
10, 11] for epidemic processes. This is identical to the SI dy¬ 
namics, except that now an infected node can revert to being 
susceptible at a (normalised) rate of ‘1’. The absorbing state 
for the SIS epidemic is when all nodes become susceptible; 
the time taken to reach this state is called the extinction 
time. We are interested in networks where this extinction 
time is exponentially large (i.e., where n is the 

number of nodes in the network). 

A key assumption in both models is that the infection rate 
for a susceptible node scales linearly in the size of the infected 
neighborhood', formally, a susceptible node j transitions to 
being infected at a rate P\Ij\, where Ij is the set of infected 
neighbors of j. This assumption makes the models concep¬ 
tually simple and tractable. However, it also means that in 
these models, high-degree nodes may get infected at unre¬ 
alistically high rates. In several real-life settings, this is un¬ 
reasonable; for example, in social networks, users broadcast 
content to all their neighbors, but attention-span constraints 
limit consumption of content to only a few neighbors. Simi¬ 
larly, in sensor networks using epidemic protocols, nodes can 
broadcast easily, but have download bandwidth constraints, 
and so can listen to only some of their neighbors. Refer to 
Section 1.2 for more examples. 

To modify epidemic dynamics to fit such settings, we need 
to assume that nodes have bounded susceptibility, i.e., a max¬ 
imum rate at which they can get infected, which is indepen¬ 
dent of network size. We can capture this via the notion 
of an infection-profile- a susceptible node j, we assume, is 
infected at a rate 4?(/3|7jj), where <!>(■) is a bounded function. 
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Now, in order to use epidemic models to understand real-life 
processes, one hopes that the behavior of these epidemics is 
somewhat insensitive to the infection-profile 4>. This moti¬ 
vates the central question we address in this paper: 

How robust is SI/SIS epidemic behavior under an 
additional bounded susceptibility constraint? 

The answer turns out to depend on the dynamics and 
network structure in a surprisingly subtle way. In particular, 
bounded susceptibility exposes a sharp dichotomy between 
two classes of networks of practical relevance: 

• Star-like networks, i.e., sparse, structured graphs, with 
a small number of edges, mostly concentrated among 
a dense core, with all other nodes at a short distance 
from the core. These are good models for engineered 
networks, which are designed to be highly connected 
while using only a small number of edges. 

• Clique-like networks, i.e., dense, uniform graphs, with 
edges spread evenly among nodes. A good example is 
the connected-regime Erdos-Renyi graph (i.e., G{n,p) 
with p = fl(logn/n)). These are good models for 
organically-formed networks. 

In prior work, the distinction between the two classes is sug¬ 
gested, but not explicitly made; we formalize this classifica¬ 
tion in Section 2.1. It turns out that this is crucial to our 
study of bounded susceptibility, as follows: 

For a fixed virulence /3, the spreading time in the SI dy¬ 
namics on star-like networks is small (O(logn), under fairly 
weak conditions on the non-core subgraph); in clique-like 
networks, the spreading time is 0(1)^. In case of the SIS 
epidemic, both clique-like networks and star-like networks 
experience exponential extinction times [8]. Thus, in a sense, 
both star-like and clique-like networks are well suited for epi¬ 
demic processes - this behavior, however, no longer holds 
under bounded susceptibility. 

1.1 An Overview of our Results 

We refer to our modified epidemic models, incorporating 
bounded susceptibility, as the Shs and SISba dynamics. As 
mentioned, we assume that in both models, a susceptible 
node j gets infected at a rate '^(/Sl/^l), where Ij are the 
infected neighbors of j; in the Sha model, infected nodes 
remain so forever, while in the SISba model, they recover to 
become susceptible at a rate 1. We assume that $(0) = 0, 
and $(■) is concave, non-decreasing, and bounded above by 

some <I>niax = ©(!)■ 

SIba Dynamics: In the Sha model, we are interested in 
finding the (expected) spreading time, under the assumption 
that both P and ‘I’max are 0(1). In this context, we show: 

• For clique-like networks, the spreading time switches 
from 0(1) under SI dynamics to f2(logn) under the 
Sha dynamics. In particular, for a G{n,p) graph in the 
connected regime (i.e., p ^ uj (logn/n)), the spreading 
time switches from o(l) to f2(logn) under the Sha 
dynamics. 

• For star-like networks, the spreading time is 0(logn) 
under both SI and Sha dynamics (under some as¬ 
sumptions on the network structure). Thus, imposing 
bounded susceptibility does not significantly alter the 
epidemic behavior. 

Thus, the SI model is highly sensitive to bounded sus¬ 
ceptibility in clique-like networks - note that the spreading 

^Special cases of this result can be found in [2, 4] 


time has switched from vanishing with network size to ac¬ 
tually scaling with network size. On the other hand, it is 
insensitive in star-like networks. 

A moment’s thought suggests that this behavior is indeed 
reasonable - we have greatly reduced the (potential) suscep¬ 
tibility of high-degree nodes, which should hinder the spread 
of the epidemic. Furthermore, since high-degree nodes are 
more affected, it seems natural that bounded susceptibility 
should adversely affect dense graphs much more than sparse 
graphs. Surprisingly, the latter intuition is completely re¬ 
versed in case of the SISba model. 

SISba Dynamics: Recall that in the SISba model, we want 
to find the (expected) extinction time under /3, 4?inax = 0(1) 
and recovery rate 1. In this context, we show the following: 

• For clique-like networks, the extinction time is expo¬ 
nential (i.e., exp(f2(poZi/(n)))) both with and without 
the bounded susceptibility assumption. In fact, this 
robustness property holds for any network which con¬ 
tains a polynomial-sized subgraph which resembles a 
random G{n,p) graph in the connected regime (i.e., 
with p ~ cu (log n/n)). 

• However for star-like networks, we now observe a sharp 
transition in extinction time. While star-like networks 
experience epidemics with exponential lifetimes under 
the SIS process [8, 9], we prove that this is not true in 
the case of SISba infections. In particular, for any net¬ 
work which can be decomposed into a small core, plus 
an outer sparse subgraph, the SISba epidemic lifetime 
is sub-exponential for any virulence P = 0(1)- 

This reversal of results may seem paradoxical - there is 
however an intuitive explanation for this. Note that the 
absorbing states in the SI and SIS epidemics are exactly the 
opposite (all-infected vs. all-susceptible). Moreover, in the 
SI model, we want the epidemic to move rapidly towards the 
all-infected state, while in the SIS, we want to deter it from 
reaching the all-susceptible state. Rapid transition to the 
all-infected state in the SI dynamics is aided by having many 
high-degree nodes - this is hampered by imposing bounded 
susceptibility. On the other hand in the SIS dynamics, high 
degree nodes impede transition to the all-susceptible state 
- having many such nodes helps distribute this effect, but 
too few high-degree nodes means that these are now critical 
bottlenecks. Under bounded susceptibility, these bottleneck 
nodes are the most affected, resulting in a sharp transition 
in the epidemic behavior. 

To make this intuition more transparent, in Section 2.2 we 
demonstrate the effect of the various dynamics in the clique 
and star networks, which are in a sense the extreme-case 
scenarios. Our main contribution in this work is to show 
how these phenomena extend to large classes of networks, 
occupying the space between these two extremes. In the 
process, we develop some structural lemmas relating the dy¬ 
namics of epidemics to network properties, which may be of 
independent interest. 

Technical Contributions. Most of our arguments involve 
standard Markov chain machinery, along with careful use of 
measure concentration. Some of the ideas we use are of fairly 
recent origin, in particular, the sharp expressions for extinc¬ 
tion time in SIS models via embedding in an ergodic Markov 
chain (Lemma 1, which we adapt from Lemma 8 in [10]), 
and characterization of the SI spreading time in clique-like 
graphs, which is similar to results for first-passage percola- 
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tion in random graphs [4]. However, unlike [4], our results 
(Theorems 1,3, 4 and 6) do not focus on a particular graph 
or generative model, but rather, relate the epidemic behav¬ 
ior to structural properties of the network. We believe these 
lemmas could be useful in studying other related settings. 
Finally, our results follow from elementary arguments, and 
hold for finite networks, with explicit error bounds in terms 
of network size. 

We note that our classihcation of networks into clique-like 
and star-like does not cover all networks. In particular, there 
are networks which are a mix of the two, and where our re¬ 
sults may not be tight; sharper techniques could perhaps be 
used to bridge this gap. Furthermore, our techniques do not 
immediately extend to other epidemic models, in particu¬ 
lar, the SIR dynamics (see [12, 13] for details) - this is an 
interesting avenue for future exploration. 

1.2 Bounded Susceptibility: Applications in 
Real-World Settings 

Before we formally state our results, we first discuss three 
important real-life settings which exhibit bounded suscepti¬ 
bility, and point out applications of our results therein. 

Social and Economic Networks. Perhaps the most com¬ 
pelling example of bounded susceptibility in network epi¬ 
demics comes from online social networks. Sites such as 
Twitter and Facebook are built on a broadcast model - users 
upload content to the website, intending it to be broadcast 
to all neighbors. However, users are known to consume con¬ 
tent from only a small set of neighbors - this is due to a 
combination of limited attention-spans and content filtering 
algorithms, and has been observed empirically in Twitter 
data [14, 15[. Epidemic models have been used to under¬ 
stand why content ‘goes viral’ on such media, however ex¬ 
isting works do not account for bounded susceptibility. 

Non-linear propagation models are also used in economics 
to model network externalities in the spread of opinions, 
technologies, etc. For example, the adoption of a new tech¬ 
nology on a network is often subject to diminishing returns 
due to wearing out of novelty with increasing number of 
recommendations [16] - this naturally suggests a concave 
infection-profile, and further, bounded susceptibility. 

Human-Disease Epidemiology. Network epidemic mod¬ 
els are widely used to study human epidemics, and also com¬ 
puter viruses. Here again, bounded susceptibility is often a 
key feature of real epidemics, due to spread of awareness 
and adoption of preventive measures. Non-linear infection 
models have been used to capture the effect of risk percep¬ 
tion in the spread of epidemics [17, 18[, by combining the 
effect of infection spread and the protection measures taken 
by the susceptible individual. However these works focus on 
approximate and mean-field treatments for the problem. 

Our work provides both intuition and technical tools for 
studying the effects of bounded susceptibility on epidemics. 
One important application of our findings is in the design 
of vaccination and quarantining policies on networks - such 
policies help reduce the susceptibility of nodes to infection, 
but usually nodes retain some susceptibility. A closely re¬ 
lated design question is of distributing antidotes in a network 
[11] - there, the epidemic is controlled by increasing the re¬ 
covery rate of a node, whereas our control is via limiting the 
susceptibility of a node. 


Distributed Computing on Networks. A third applica¬ 
tion of epidemic processes is as a primitive for designing 
low-complexity distributed algorithms. The SI dynamics is 
the basis of the flooding algorithm, which is used for broad¬ 
cast messaging in networks, distributed database synchro¬ 
nization, etc [19]. The SIS dynamics is sometimes used to 
ensure persistent storage in distributed storage/sensing ap¬ 
plications with disk failures [20]. Existing work focuses on 
using the standard SI/SIS dynamics, or modifies them to 
impose bounded influence, i.e., a bound on the upload band¬ 
width of a node. 

However many of the above settings, in reality, exhibit 
bounded susceptibility. Modern distributed/P2P systems 
are based on a virtual network architecture, where nodes 
upload content to a central location, from where it is down¬ 
loaded by its neighbors. Thus, broadcasting is easy, but 
downloading is bandwidth constrained. This is also true 
for sensor networks - a node needs to broadcast a message 
once, but it is limited in how many neighbors it can listen 
to, due to a fixed download bandwidth. Thus, bounded sus¬ 
ceptibility is natural in these settings. We discuss all these 
applications in Section 2.4, after presenting our results. 

1.3 Related Work 

Epidemic processes on networks have been studied across 
many disciplines; readers interested in more details regard¬ 
ing epidemics are referred to several excellent books on the 
subject [1, 12, 21, 22]. Specifically, there is a vast literature 
to characterize spreading time in various contexts for SI pro¬ 
cesses [1, 2, 3, 4], and phase transitions/extinction time for 
SIS processes [6, 7, 8, 9, 10, 11]. Phase transitions for SIR 
processes are available in [13]. 

Non-linear epidemic models have been explored in the past 
in various contexts. A large body of work is devoted to set¬ 
tings with bounded influence, where the total rate at which 
a node infects its neighbors is bounded [23, 24], but the 
infection-profile '!>(•) is still linear. One popular model with 
non-linear '!>(•) is the Bootstrap Percolation, and related 
variants [25, 26]. Several authors have also studied non¬ 
linear infection models at the population level using mean 
field approximations [17, 18] - these works show the exis¬ 
tence of epidemic thresholds under different non-linear mod¬ 
els. However, to the best of our knowledge, ours is the first 
work that explicitly characterizes epidemic behavior under 
bounded susceptibility. 

2. MAIN RESULTS AND DISCUSSION 

We now present and discuss our main results. Our gen¬ 
eral theorems are presented in Section 2.3, followed by a 
discussion of their application in Section 2.4; before that, in 
Section 2.2, we specialize our results for the clique and star 
networks, and give an outline of the proofs for these cases. 
The complete proofs are somewhat technical - we outline 
the main proof ideas after each theorem, deferring complete 
proofs to Section 3. 

2.1 System Model 

The SI/SIS Dynamics 

The SI, SIS dynamics models we use are standard in litera¬ 
ture [8, 10]. We consider a graph G{V, E) with n nodes (ver¬ 
tices). The various processes [SI, Shs, SIS, SISbs) evolve 
on this graph in continuous time. Associated with each node 
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i € {1, 2,... ,n} is a random process Xi{t) G {0,1}. Here, 
‘0’ corresponds to the susceptible state and ‘1’ corresponds 
to the infected state (the entire state vector is denoted by 
X(t)). For the SI and SIS models, nodes in state ‘1’ (in¬ 
fected nodes) attempt to infect a neighbor independently at 
a rate /?, i.e., after a i.i.d random interval drawn from an 
Exponential{l3) distribution. In the SI dynamics, an in¬ 
fected node never recovers. In contrast in the SIS dynam¬ 
ics, an infected node returns back to the susceptible state at 
rate 1 (i.e., after an Exponential{!) time). 

We assume that the infection starts off at a single node, 
arbitrarily chosen. This is natural for the SI epidemic since 
we want to study spreading-time. For the SIS epidemic, 
our results generalize for any constant-sized initial infected 
set. For the SI dynamics, the only absorbing state is the 
all-infected state (i.e., Xi = 1 Vi), and the spreading-time 
is defined as Ts = inft{Xi{t) = 1 Vi}. In contrast, for the 
SIS dynamics, the only absorbing state is the all-susceptible 
state (i.e., Xi — 0 Vi), and the extinction-time is defined as 
Te = inft{Xi(t) — 0 Vi}. We also assume throughout that 
/? is a constant, independent of n; however, our analysis 
extends for settings where /3 is a function of n. 

Infection-Profiles and Bounded Susceptibility 

A critical assumption in both SI and SIS models is that 
the infection rate for a susceptible node scales in proportion 
to the number of infected neighbors. In order to study the 
effect of non-linearity in the infection-rate, we introduce the 
notion of an infection-profile: we assume that a susceptible 
node j is infected at a rate 'I‘(/3|7j|), where Ij is the number 
of infected neighbors of node j (at some time f; we suppress 
the dependence on t for ease of notation). Thus, in both the 
SIba and the SISba dynamics, we have: 

Xi : 0 —>■ 1 at rate 'h I /? E > 

V (ij)eB J 


We can also define an alternate notion of isoperimetry based 
on vertex expansion, as follows: 

where T{S) is the neighborhood of S, i.e., the set of nodes in 
S^ with at least one neighbor in S. It is easy to check that 
for any fixed m, we have rj^{G) > rjf^iG). 

The isoperimetric constants characterize bottleneck sets, 
i.e., those containing the smallest fraction of the potential 
edges in the cut. For a clique, it is easy to check that r]{G) = 
0(n), which is the highest possible. To show that epidemic 
properties hold over a large class of graphs, we want to admit 
graphs with smaller isoperimetric constants. This motivates 
the following characterizations for clique-like networks: 

Definition! (Clique-like Networks). We define a 
graph G to be clique-like if it satisfies one of the following: 

(A) IfrjifG) = H(logn). 

(B) If for some m = n'^, a £ (0,1), we have rfflG) = tu(l). 

The two definitions arise from analyzing the two epidemic 
dynamics on graphs. Although not equivalent, they both 
characterize networks with good expansion properties. To 
demonstrate their wide applicability, we show that random 
(Erdos-Renyi) graphs in the connected regime (i.e., G{n,p) 
with p ^ Q. (logn/n)) are clique-like under both definitions. 

On the other hand, star-like networks intuitively are those 
which can be partitioned into a dense core and a sparse 
periphery. Here, the natural notions of density turn out 
to be the average degree davglG), and the spectral radius 
p{G), i.e., the largest eigenvalue of the adjacency matrix, 
which is closely related to the node degrees (in particular, 
davg{G) < p{G) < dmafiG). Formally, we write V = VcUVp, 
where 14 denotes the core and Vp the periphery. As for 
clique-like networks, we have two characterizations: 


and in the SI Sts, we additionally have Xi : 1 —>■ 0 at rate 1. 

$ : R+ —>■ ]R+ represents the overall effect that the infected 
neighbors together have on a susceptible node. We assume 
that 'I?(0) = 0, and that <!>(•) is non-decreasing, concave and 
bounded. Note that these assumptions imply that := 

supj,gR 3>(a;) exists, and $(/3) > OV/3 > 0. 


Clique-Like and Star-Like Networks 

We now formalize the two classes of networks, which we in¬ 
troduced in Section 1. Intuitively, clique-like networks are 
those which are dense (i.e., where |i5| = oj(n)) and have 
edges evenly distributed across the cuts of the network. One 
way to formalize this is via the expansion properties of the 
underlying graph. There are different ways to quantify the 
expansion of a graph. For example, the edge expansion of a 
graph is characterized by the generalized isoperimetric con¬ 
stant: 

f?m(G) := inf 0<m<Ln/2j, 

5CV,|S|<m |o| 

where, for any set A £V, S(A) = {(*, j) ^ E s.t. i £ A,j ^ 
A\ are the edges in the cut defined by A. The case of m = 
n/2 is referred to as the isoperimetric/Cheeger constant: 


V{G) := 


inf 

SCU.|S|<n/2 |y| 


Definition 2 (Star-like Networks). We define a 
graph G to be star-like if it satisfies one of the following: 

(A) If \Vp\ = Q{n) and the average degree of nodes in 
the periphery is 0(1). In addition, the diameter of the 
graph is O(logn). 

(B) //|I41 = O{polylogn) and the subgraph G{Vp) induced 
by the periphery has spectral radius p{G{Vp)) = o(l). 

Again, the class of graphs that are star-like by the two def¬ 
initions above are not equivalent but broadly include net¬ 
works which have a sparse periphery and a small (possibly 
dense) core. We derive our results for SI dynamics based 
on definition [A) of clique-Iike and star-like networks, while 
the results for SIS dynamics are based on definition (B) for 
both the networks. 

2.2 Two Examples: The Clique and the Star 

Before presenting our full results, we first illustrate them 
by considering our two representative graphs -- the clique 
and the star. This comparison brings out the dichotomy 
between the two network classes. Further, it allows us to 
build intuition for our main results, since these special cases 
admit simple proofs. In all these results, we use standard 
notation (O, 0, H , o and lv) to characterize the scaling be¬ 
havior with the graph size, n. We also use the notation 
H = {1,2,... ,n}, and Hn = f 
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SI - Shs Dynamics 

Consider first the SI/Sits epidemic on K^, the clique on n 
nodes. We now show that the spreading-time is vanishingly 
small under the SI dynamics, while it scales with n under 
the SIba dynamics. 

Proposition 1. For the SI/Shs epidemic on start¬ 
ing at an arbitrary node, 

1. Under the SI dynamics, E[ts] = 0( 

2. Under the SIbs dynamics E[ts] = 0(logn). 

Proof. The first part is known from previous results [4] 
and follows from standard Markov-chain arguments. By de¬ 
noting N{t) to be the number of infected nodes at time t, and 
analyzing the steady-state of the resulting (one-dimensional) 
CTMC, we immediately get the result. 

Next, for the bounded-susceptibility spreading process, 
SIbs, we again consider the process {N{t)}, i.e., the total 
number of infected nodes. {N{t)} has transition rates as 
follows: 

z z -I 1 at rate 4>(/32)(n — z) Vl< 2 <n — 1, 
and hence the expected spreading-time is given by: 

I I Hn-l 

^ 'h(/3fc)(n - fc) “ ^ (n - fc)$mQa: ^max’ 

where the inequality follows since ^(pk) < $rnax- From the 
same equation, we also have: 

ip-r 1 _ ^ _ I _ ^ Hn-l 

^^HPk){n-k) - ^p)’ 

since $ is a non-decreasing function. Finally, since both 
4>(/3) and 4>inax are 0(1), and using In(n-l-l) < Hn < ln(n)-|- 
1, we get the second assertion. □ 
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Figure 1: Simulated epidemic spreading-time in a 
clique under SI and SIbs, averaged over 500 runs for 
P = Q.l and ^max = 0.5. Note that E[rs] switches from 
o(l) under the SI dynamics, to f2(logn) under SIbs 
(decreasing to increasing on log — log plot). 



SIS - SISbs Dynamics 

Next, we consider the SIS/SISbs epidemic. Here, it is 
known [8] that the expected extinction-time is exponential 
in cliques - we will prove that it is exponential even for the 
SISbs dynamics. 

We first state and prove a lemma that gives a closed-form 
expression for expected time to absorption in a birth-death 
Markov chain. The lemma is adapted from [10, Lemma 8] 
- it re-obtains the result for SIS epidemic on cliques, and 
also extends to our new result for the SISbs epidemic. 

Consider a birth-death CTMC {U{t)} on {0,1,..., n}, 
with the following transition rates: 


Thus there is a sharp change in the spreading-time be¬ 
tween the SI and SIbs epidemic processes in cliques. This 
is clearly visible in Figure 1, where we simulate and plot the 
spreading-time under the SI and Shs dynamics on a clique, 
as a function of network size n. Further, in Theorem 8 in 
the Appendix, we derive concentration bounds showing that 
the above result holds with high probability. 

On the other hand, the same behavior is not seen in a star 
network, as we show via our next result: 

Proposition 2. In a star network with a single hub and 
n—1 leaves, under both the SI and SIbs epidemics, we have: 

E[ts] = 0(logn) 

Proof. Consider first the SI process on a star network. 
We denote the hub as node 1 and define: 

ri = inf{Node 1 is infected}, 

i.e., the first time that the hub gets infected. Since the epi¬ 
demic starts at a single (arbitrary) node, we can split the 
spreading-time as: Ts = ti + (ts — ti). Next, we have that 
E[ti] = 0(1) and Ts — ri corresponds to the maximum of 
m G {n — 2, n — 1} i.i.d random variables, each drawn from 
a Exponential{P) distribution; thus E[rs — n] = 2^ r.., 
0(logn). Finally, note that no susceptible node ever has 
more than 1 infected neighbor - hence the Shs model is 
identical to the SI model with P' — ^{P) = 0(1)- There¬ 
fore, the result also holds for SIbs epidemics. □ 


i ^ i + 1 at rate qi,i+i Vl<i<n — 1, 
i —>■ i — 1 at rate V1 < i < n. 

Note that 0 is the only absorbing state for the above process. 
Let T{i,j) := inf{t : U{t) = j\\U{0) = i}. Then we have: 


Lemma 1 (Lemma 8 from [10]). 


E[r(i,0)] 


1 

qi,o 


i+En 


Qi,i+l 


For the sake of completeness, we provide a proof of this result 
in the Appendix. We now use this to get lower bounds on 
the expected extinction-time in cliques. 


Proposition 3. For the SIS/SISba epidemic on clique 
Kn starting at an arbitrary node, under both the SIS and 
SISbs dynamics, log(E[Te]) = fl(n). 

Proof. We prove the statement for the SISbs dynamics; 
the claim for the SIS dynamics follows similarly (and is 
also known from previous work [8]). As in Proposition 1, it 
is sufficient to consider the total number of infected nodes 
N{t) at time t. {A(t)} has transition rates as follows: 

z —>■ z -h 1 at rate $(Pz)(n — z) Vl< 2 <n — 1, 
z ^ z — 1 at rate 2 V1 < 2 < n. 

Let ko = inf{fc G N : ^(Pk) > ^■max/’i/- Since "F is concave 
and $(0) = 0, we have ${pk) > -^^(pko) > Vfc < 
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ko. Note that fco = 0(1)- Now, from Lemma 1, we have: 


E[re] = 1 + n 


k = 2 i=l 

kn k — 1 


^1+En 


k=2i=l 

n fcg — 1 


— i) 

i + 1 

^max iijl i') 
2/1)0 Z “h 1 


^max ^(tT, — i) -p-r ^rnax ('^ “ j) 
2^ LL O.hn i + 1 11 


= 1+E 

k = 2 
fco! 


2fco j + 1 
fc = fco + l i=l j = kQ 

(fc-l)! /n\ 


i + 1 


2 fco 


+ ■ 


nk: 


E 


> 


0 A:=fco + l 
2fco! 


r7<t) 

'^^max^Q k = l 


■E 


n\ 

fc 

n\ / 4?^ 


2fco! 




1 + 


$r, 


-1 . 


Since 4>maa;,fco are 0(1), log(E[re]) = Q,{n). □ 


The more interesting phenomenon in the SIS model oc¬ 
curs in the star network. As in the clique,it is known [ 8 , 9] 
that the expected extinction-time in star networks is expo¬ 
nential. However, nnder the SISba dynamics, we show that 
it is sub-exponential. More precisely, we show that the epi¬ 
demic survives in a star network for a time that is at most 
polynomial in n. 


Proposition 4. For an SISbs epidemic originating at 
any node in a star graph, E[re] = 

This change from exponential to polynomial infection time 
is clearly seen in simulation results shown in Figure 2. The 
proof of this proposition is somewhat technical, so we pro¬ 
vide only a proof sketch here. The high-level proof structure 
is similar to existing analysis of the SIS epidemic on a star 
[ 8 , 9] - the novelty lies in showing that bounded susceptibil¬ 
ity causes a sharp transition in the epidemic lifetime. 

Proof Outline. We consider cycles of the epidemic evo¬ 
lution - each cyclical epoch starts off with the hub infected, 
then recovering to become susceptible, and finally getting 
reinfected, thereby starting the next cycle. The time from 
the start of an epoch to when the hub first becomes sus¬ 
ceptible is an i.i.d Exponential{l) variable - at this time, 
some subset of leaf nodes are infected. Subsequently, ei¬ 
ther the hub is re-infected (i.e., a new epoch starts), or all 
these leaf nodes recover before the hub is re-infected -(i.e., 
the epidemic becomes extinct). The crux of the proof lies in 
estimating the probability that, under the SISba model, the 
leaf nodes recover before the hub is reinfected. Given this, we 
then show that the total number of such epochs is stochas¬ 
tically dominated by a geometric random variable, whose 
mean is polynomial in n - this gives the above result. □ 

In Section 2.3, we state a much more general version (The¬ 
orem 6 ) of this result, extending it to star-like networks. The 
complete proof of both results are given in Section 3. 


2.3 Main Results 

We now state our main results, which generalize results 
from the previous section to clique-like and star-like graphs. 



Figure 2: Simulated epidemic extinction-time in a 
star under SIS and SISba, averaged over 500 simula¬ 
tion runs for /3 = 0.1 and ^max = 0.5. Note thatE[re] 
switches from exponential under SIS to polynomial 
under SISba (polynomial to linear on log —log plot). 


2.3.1 SI — SIba Dynamics 

Recall that, for a graph G, the Cheeger constant is given 
by ri{G) — infscv,|S|<n /2 We now have the following 

theorem, which allows us to bound the spreading-time of the 
SI epidemic in clique-like graphs: 

Theorem 1. Given underlying graph G with Cheeger con¬ 
stant rj{G), the spreading-time for an SI epidemic with vir¬ 
ulence jl originating at a single arbitrary node, obeys: 

1 / (2 + o(l))log(n) 

--WG)-' 

The proof of this theorem, presented in Section 3, is sim¬ 
ilar to that of Proposition 1. Note that for the clique, we 
have 17 (G) = n/2. Plugging this into the above bound im¬ 
mediately gives back Proposition 1. However the theorem 
captures in finer detail the impact of expansion (as captured 
by giG)) on spreading-time. As a corollary to the theorem, 
we have the following result for clique-like networks: 

Corollary 1. In a network that is clique-like as per Defin- 
tion 1(A) (i.e., g{G) = H(logn)^, the expected spreading 
time of an SI epidemic is 0(1). 

We note also that although Theorem 1 bounds the expected 
spreading-time, we can extend it to get a concentration 
bound for the spreading-time; refer to Theorem 7 in the 
Appendix for details. To see how the result can be used 
to characterize the spreading-time in a variety of clique-like 
graphs, we present the following corollary for the case of 
dense random graphs, i.e., Erdos-Renyi graphs in the con¬ 
nected regime: 

Corollary 2. Let G be an Erdds-Renyi graph G(n,p), 
with p > ” . Then with probability at least 1 — , we 

have r)(G) > Therefore, for an SI epidemic, 

< (8 + o(mosn 
' - I3np 

with probability at least 1 — . 
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Note that by the Borel-Cantelli Lemma, the above bound on 
the expected spreading-time holds almost surely as n —^ oo. 
The proof is provided in Section 3. 

Putting the pieces together, we can now compare the SI 
and SIbs processes in dense random graphs: 

Theorem 2. Let G be an Erdds-Renyi graph G{n,p), with 
p > ^ . For the SI/SRs epidemic on G, starting at an 

arbitrary node, we have: 

1. Under the SI dynamics, E[rs] = O with proba¬ 

bility at least 1 — 

2. Under the SIbs dynamics E[ts] = fl(logn). 

Proof. The first claim follows from Corollary 2 and The¬ 
orem 1. For the second, note that given any graph G, and G' 
obtained by adding edges to G, Ts{G') <st rs{G) - this can 
be shown via a standard coupling argument. Thus, for any 
graph G, we can always lower bound the spreading-time un¬ 
der the SIbs dynamics by that on the clique K^. The second 
claim thus follows from Proposition 1. □ 

Next we turn to the case of star-like graphs. To charac¬ 
terize SI epidemics on such networks, we first state a lemma 
that gives a lower bound on E[ts] in terms of average degree: 

Lemma 2. Consider a graph G with \ V\ = n, and an SI 
epidemic on G, starting at some node vi. Let A I- V be 
a subset of nodes containing vi, and let be the average 
degree of the nodes in V \ A, i.e., dA<= = Then: 

E[r.] > 1^1) 

PdA- 

For example, for the star graph with vi chosen to be the hub, 
we have d^^^yc = 1, and so we get back the lower bound 
in Proposition 2. We now use the above lemma to show 
that imposing bounded susceptibility does not significantly 
change the spreading-time in star-like networks: 

Theorem 3. Consider a network G of diameter O(logn) 
that is star-like by Definition 2{A). In other words, there 
exists a partition of the nodes into a core 14 and periphery 
Vp, parametrized by: 

• Core size: \Vc\ = m is at most a constant fraction of 
nodes (i.e. m < cn for some c < 1). In other words, 
| 14 |=n(n). 

• Sparse periphery: dp — 0(1), i.e., the average degree 
of nodes in the periphery is constant (or lesser). 

Then both under the SI and SIbs dynamics, E[rs] = 0 (logn). 

Note that the sparsity constraint is on the average degree 
of nodes in Vp, not in the subgraph induced by Vp, i.e., it 
also accounts for the edges in the cut (5(Vp). 

Proof. First consider the SI dynamics, starting at some 
arbitrary node ui. We want to derive a lower bound on 
the spreading-time using Lemma 2 - to do so, we choose 
A = 14U{ui}, and hence A‘^ = V),\{ui}. Now using Lemma 
2 and the fact that 1141 ~ poly{n), and that both /3 and dp 
are 0(1) we have that E[ts] = f2(logn). 

To obtain an upper bound on the spreading-time, we under¬ 
dominate the spread by considering epidemic spread on the 
shortest-path spanning tree of the network rooted at the 
initial infected node. The resulting depth is at most the 


diameter of the tree. Thus, the spreading-time is stochas¬ 
tically upper bounded by the maximum of n i.i.d random 
variables distributed as sum of diam{G) number of exponen¬ 
tial random variables with rate p. Finally, using standard 
concentration results for the sum of exponentials (refer [3]), 
we have that E[ts] = O(logn). 

For SIbs epidemic, we can first over-dominate the epi¬ 
demic spread by converting the network into a clique - we 
have from Proposition 1 that the time taken to infect all 
nodes is fl(logn). On the other hand, we can again under¬ 
dominate the spread by considering an Shs process on the 
minimum spanning tree which is equivalent to an 57 process 
of rate $(/3). Consequently, E[rs] = 0(logn). □ 

2.3.2 SIS — SISbs Dynamics 

Next we look at the SIS — SISbs dynamics on clique¬ 
like and star-like graphs. Recall that in this case, we are 
interested in the extinction-time Te. In particular, we are 
interested in characterizing networks in which the epidemic 
is persistent, i.e. the extinction-time is exponential in n (i.e., 
E[re] ~ n(eP°'^("))). 

Existing works [8, 9, 10] demonstrate that the virulence 
thresholds for short-lived (E[re] ~ O(logn)) SIS epidemics 
depend on the spectral radius, p(G), while thresholds for 
persistent epidemics are related to the generalized isoperi- 
metric constants. Though the results for the SIS epidemic 
depend on edge-expansion isoperimetric constant 774 (G), the 
vertex-expansion isoperimetric constant rff, (G) turns out to 
be the appropriate notion for the SISbs epidemic, as follows: 

Theorem 4. Consider an SISbs epidemic on graph G, 
with virulence P and infection-profile $(•). Suppose form = 
n“ for some 0 < a < 1, we have that: 

HPWmiG) > 1 , 

then E[Te] = \ i.e., the SISbs epidemic is persistent . 

Proof Outline. The proof is similar to that of Propo¬ 
sition 3. In particular, we construct an under-dominating 
CTMC for the number of infected nodes under SISbs. Next, 
we use Lemma 1 to get an expression for the extinction¬ 
time, and show that the bound is indeed exponential in n if 
^{P)Vm{G) > 1. For complete proof, see Section 3. □ 

An immediate consequence of this result is the following 
corollary for clique-like networks. 

Corollary 3. An SISbs epidemic is persistent in any 
network belonging to the class of clique-like networks as per 
Definition 1{B). 

Returning to the SIS epidemic, the strongest known condi¬ 
tion for a persistent SIS epidemic is Pr](^{G) > 1 [8]. Now 
since rj(r,{G) > r](n{G), this implies that Prf(^{G) > 1 is 
also sufficient for a persistent SIS epidemic, i.e., without 
bounded susceptibility. Theorem 4 shows that the condition 
for a persistent SISbs epidemic is closely related. Further¬ 
more, as we show subsequently, the vertex-expansion condi¬ 
tion is strong enough to recover all the results from [8], in 
particular for the G(n,p) and power-law graphs. 

We now use Theorem 4 to demonstrate that both SIS and 
SISbs epidemics have exponential extinction-times in dense 
random graphs: 
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Theorem 5. For both the SIS and SISf,s epidemics on 
a G(n,p) graph withp > g, large enough n, with 

probability at least 1 — uie have, 

log (E[re]) = r 2 (n). 

Proof Outline. We show under the conditions of the 
theorem that ^{/3)ril^{G) > 1 with probability > 1 — if 
m = n°‘ for any constant a G (0,1). The resnlt then follows 
using Theorem 4. Complete proof is given in Section 3. □ 

Note also that any graph that contains a large clique-like 
subgraph is prone to an exponentially long lasting infection, 
even under bounded susceptibility, if the initial infected node 
belongs to the subgraph. Clearly a clique, or a realization 
of an Erdos-Renyi graph in the connected regime experi¬ 
ences long lasting infection - however such dense random 
subgraphs appear in other graph families including ‘power- 
law’ or scale-free networks. We discuss this further in Sec¬ 
tion 2.4 (in particular. Corollary 5). 

More generally, these results show that large dense net¬ 
works experience persistent epidemics even under bounded 
susceptibility due to the presence of large cligue-like sub¬ 
graphs. This is in sharp contrast to existing work (for exam¬ 
ple, [9]) which suggests that persistent epidemics in such net¬ 
works arise due to the presence of high-degree nodes (‘hubs’). 

In fact, it turns out that the presence of such ‘hubs’ do not 
help - we took a step towards showing this in Proposition 
4, where we showed that a star graph does not experience 
persistent SlSbs epidemics. We now show that this is true 
for a much larger class of star-like graphs. As before, the 
node-set V is decomposed into core 14 and periphery Vp. 
We characterize the sparsity of the periphery via the spec¬ 
tral radius p{G{Vp)) of the periphery subgraph - recall that 
davg{G{Vp)) < p{G{Vp)) < dina.x(G(Vp)). Now we have: 

Theorem 6 . Consider a network G, with the nodes par¬ 
titioned into a core 14 and periphery Vp characterized by: 

• Core size: |141 = rn. 

• Periphery sparsity: Parametrized by p{G{Vp)). 
Suppose G{Vp) satisfies ^{P)p{G{Vp)) = 1 — e, for some 
£ > 0. Then for an SISba epidemic on G,: 

log (E[Te]) = 0 (mlogn). 

Proof Outline. The high-level proof structure is simi¬ 
lar to the outline we gave for Proposition 2. We again con¬ 
sider cyclic epochs. An epoch starts with some node in the 
core infected, transitioning to all core nodes becoming sus¬ 
ceptible, and finally ends when some node in the core gets 
reinfected beginning the next epoch. Next, we bound under 
the SISbs model the probability that, in any epoch, the epi¬ 
demic dies out in all nodes of the periphery before any node 
in the core is reinfected. Finally, via stochastic domination 
by a geometric random variable, we get the above result. 
The complete proof is given in Section 3. □ 

Though we state the above theorem for a general core size, it 
is most significant in settings the core is small. For example, 
if the core is of constant size, then E[re] scales only polynomi- 
ally with n - note that this includes the star network (where 
pifGiVp)) = 0, as the periphery is fully disconnected if we re¬ 
move the hub) that we characterized in Proposition 4. More 
generally, the following result for star-like networks holds as 
a result of Theorem 6 . 


Corollary 4. The extinction time of an SISbs epidemic 
is sub-exponential in any network that is star-like as per Def¬ 
inition 2{B). 

We note that we can also derive concentration bounds, show¬ 
ing that the bound on the order of the extinction-time holds 
with high probability (see Theorem 9 in the Appendix). 

2.4 Discussion 

We now show how our results apply to the settings we 
discuss in Section 1.2. 

Social and Economic Networks. Epidemic models have 
been widely used to model the viral spread of content in 
social networks. In particular, the SIS model can be used 
to study the spread of ‘memes’, i.e., related pieces of con¬ 
tent that re-appear cyclically. Furthermore, the exponential 
extinction-time of an epidemic is in a sense connected to 
long-lasting nature of popular memes. 

However, the existence of popular memes is often regarded 
to be as a result of highly influential nodes, i.e., those with 
high degree. Such high-degree nodes are a distinctive fea¬ 
ture of power-law graphs, which are often used to model 
such networks. Our results however suggest that persistent 
epidemics arise due to large random subgraphs rather than 
high-degree nodes. We now formalize this notion by apply¬ 
ing Theorem 5 to power-law graphs. 

We consider the case of power law graphs with exponent 
2 < 7 < 3. For random power law graphs generated under 
the Chung-Lu model [27] with exponent 2 < 7 < 3 , it is 
known that they contain a n”’’ connected random graph [ 8 ], 
and further, that they contain an clique [28] for some 
ei ,£2 > 0. The SIS epidemic is known to have exponential 
lifetime in such graphs [ 8 ]; using Theorem 5, we recover this 
behavior under bounded susceptibility: 

Corollary 5. For networks generated via the Chung-Lu 
model [27], if the power-law exponent 7 G (2, 3) and ^max > 
1 + u for some u > 0, the SISbs epidemic has exponentially 
long lifetime. 

Thus scale-free networks do experience persistent epidemics 
even under bounded susceptibility, due to the presence of 
large clique-like subgraphs. This guestions the belief that 
highly-connected nodes are most responsible for viral content. 

Epidemiology. Epidemic models help understand the dy¬ 
namics of infectious diseases amongst human/animal pop¬ 
ulations, and the spread of computer viruses and worms. 
These models also help in the design of defense mechanisms 
against such epidemics - for example, policies for spread¬ 
ing antidotes, which increase the recovery rate in the SIS 
epidemic ]11]. In this regard, our work helps in developing 
a theory for partial quarantine policies, wherein adopting 
preventative measures helps reduce (and bound) the sus¬ 
ceptibility of a node. We make a distinction here between 
partial and perfect quarantine policies, wherein a node can 
be completely disconnected from the network. Achieving a 
perfect quarantine is often unfeasible and there is always 
some residual susceptibility. 

In this regard, our results paint a somewhat bleak picture 
- enforcing quarantine policies in dense networks has an in¬ 
significant effect on the epidemic lifetime. More specifically, 
as long as there is some 0 ( 1 ) residual susceptibility, the epi¬ 
demic lifetime remains exponential in clique-like networks. 
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Note also that a similar effect is seen in [11] (under linear 
susceptibility), where an effective policy needs to give an 
amount of antidote to a node which is proportional to its 
degree - this is clearly more difficult in dense graphs. 

The story is somewhat more positive in star-like networks, 
where Theorem 6 shows that in a network with a small core 
and a sparse outer network, partial quarantining affects a 
large change in the extinction-time, reducing it to polyin). 
Star-like networks are relevant for studying the spread of an 
infectious disease through human transportation networks, 
which comprise a core of long distance air-links, overlaid on a 
slower road network. In case of an epidemic, a natural policy 
is to focus on safeguarding the core cities (via quarantining). 
For this, we have the following corollary: 

Corollary 6. Consider a network which can be decom¬ 
posed into: 

• A dense core of size 0(1), which can be fully connected. 

• A sparse periphery, with p{G{Vp)) — 0(1) (no restric¬ 
tion on edges between periphery and core). 

Now suppose only core nodes have bounded susceptibility 
(via quarantining), and the unprotected periphery satisfies 
l3p{G{Vp)) < 1 — e for some e > 0. Then the epidemic is not 
persistent, i.e., E[re] = 0{poly{n)). 

Note that given fl{n) edges between the hub and the periph¬ 
ery, the SIS epidemic has an exponential lifetime. Thus, 
safeguarding hub cities does help reduce the epidemic life¬ 
time dramatically. However, distributing antidotes may be a 
much more effective policy than quarantining as it can effect 
much faster (upto 0{\ogn)) extinction-times. 

Distributed Computing Applications. Another important 
application of SI epidemics is as a primitive for algorithms 
like flooding [19], i.e., forwarding a message to all network 
nodes. In order to minimize the flooding time, one approach 
would be to make the underlying network clique-like - how¬ 
ever our results indicate that under bounded susceptibility, 
the flooding time is fl(logn), which is identical to that in a 
star-like network (for example, a hierarchical network, with 
a constant-sized core and all other nodes connected to the 
core). On the other hand, our results also show that star- 
like networks are much easier to protect from viruses which 
spread like an SIS process. Considering the tradeoff between 
performance and robustness to infection, star-like architec¬ 
tures appear more appropriate than clique-like architectures 
for distributed computing applications. 

In case of networks using SIS epidemics to ensure persis¬ 
tent data storage [20], the story is strikingly different from 
the above. Here, the presence of bounded susceptibility re¬ 
sults in a sharp difference between epidemic lifetime in star- 
like and clique-like networks. However, our results also in¬ 
dicate that it is not necessary to have very high edge densi¬ 
ties, or very well-designed networks in order to achieve long 
lifetimes - on the contrary, as long as there exists a large- 
enough connected random subgraph, the epidemic lifetime is 
exponential even under bounded susceptibility. 

3. PROOFS 

We now give detailed proofs for all our results. We use 
the standard stochastic dominance notation for 

relations between random variables. Specifically, {X <st 

Y) V2,P[A > 2] < P]y > 2]. 


Proofs from Section 2.3.1. 

Proof of Theorem 1. Let N{t) be the number of in¬ 
fected nodes at time t and Tk be the hrst time at least k 
nodes are infected, i.e., Tk — inf{t : N{t) = k}. Also, 
V fc G ]n], let Vk denote the set of all subsets of V of size k, 
and let Sk G Vk be the set of infected nodes at Tk. 

Since Ti — 0 and Tn ~ Ts, we can telescope and write 
Tg = ffkZfiTk+i — Tk). Recall that in case of the clique in 
Proposition 1, we had that E[rfc+i — Tk] = pk{n-k) ■ This is 
no longer true, but we can upper bound the time in terms 
of the isoperimetric constant piG), as follows. 


E]Ti;+i — Tk] = E]ri;+i — Tk]Sk = Sfc]P]S';c = Sfe] 

1 


= E 


L/5|5(5fc)|J ■ 

By the definition of p(G), 

Jp(G)k if fc < n/2 

^r;(G)(n — k) if k > n/2. 


\S{Sk)] > 

Therefore, we have 

Ln/2J 


n — 1 


E]ts] < 


Pp{G)k 


+ E 


1 


k=l : [ra/2J+l 

1 ^/ 2 J 


Pv{G){n-k) 


<2E 


k=l 


pp{G)k 


< 


(2-b o(l)) log(n/2) 

Pn{G) 


This completes the proof. □ 


Next we show that a G{n,p) in the connected regime has 
a large isoperimetric constant with high probability. 

Proof of Corollary 2. Let graph G be a realization 
of a G{n,p) random graph with p > 2i22SJl^ YIe claim that 

P [ 77 (G) <np/4] < -^. 

n‘‘ 

In other words, for all but a vanishing fraction of graphs (of 
measure < f/rd), we have the required bound on 77 (G). To 
show this, we write 


P [ 77 (G) < np/4] 


u u ]5(sfc)] < nfcp/4 

fe<[n/2j SfceVfc 


< 


u u ]5(sfc)] < k{n - k)p/2 

fc<[n/2j SfcGVfc 


- E E 

k< [n/2J SfcSVfc 


We can now use the following Chernoff bound - for any r.v. 
B taking values in [0,1], we have 

P [B < (1 - 5)E[B]] < exp . 


Since ](5(sfe)] ~ Binom{k{n — k),p), 


Ln/2J 

P [ 77 (G) < np/4] < Y 

k = l 


exp 


-1 


k{r 



< n 
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Using the condition p > ” , we get our result. □ 

Next, we turn to Lemma 2, wherein we derive a lower bound 
on the spreading-time in terms of the average degree of a 
graph. First, we need the following lemma: 

Lemma 3. Let {Xi}^i be independent exponential ran¬ 
dom variables with rate Oi (i.e., Exponential{ai)) such that 

be i.i.d exponential random 
variables with rate Ajm. If Xmax = max{Xi, X 2 ,..., Xm} 
and Lmax = max-[Vi, L 2 ,. •., Xm}, then Xmax ^st L^ax. 

Proof. For stochastic dominance, we want to show that 
P [Xmax < *] < P [Umax <x\'ix > 0. For any X > 0, de¬ 
fine Qx(a.) — “ exp(—ai®)). We now show that a* = 

[Ajm, Ajm,... ^Ajm) maximizes over the set ^ = {a : 

For rn = 2, one can check that a* = 
(j 4/2, A/2) is the unique maximum. Thus, Qa,(a*) > Qx(a) 
for any a 7 ^ a*. If m > 2, suppose that a G A and a yf a*, 
maximizes Qx over A. Then 3i,j such that hi yf ay. Then, 
for a equal to a except with the i,j elements replaced with 
{hi -f ay)/2, we have Qx{si) > Qx{sl), which is a contra¬ 
diction. Therefore, a* = {A/m, Afm,..., A/m) maximizes 
Qx over the set A, and P [Umax < a;] = Qx{a*) > Qx{a) = 
P [Xmax < x] for any a; > 0 . □ 


Using this, we now prove Lemma 2. 


Proof of Lemma 2. Let the nodes be numbered from 1 
to n such that A = { 1 , 2 ,...,|A|} and let d = (di, d 2 , ■ ■ ■ ,dn) 
be the degree sequence of G. At any time, the rate at 
which a susceptible node i is infected is upper bounded by 
fddi. Therefore, Ts >st max{X|A|+i,..., X„} where {Xi} 
are independent exponential random variables with mean 
l/Pdi. By Lemma 3, Ts >st Y, where Y is the maximum 
of n — |A| i.i.d. exponential random variables with mean 

yi—|A _ 1 

L)"=|A| + 1 


E [ra] > E [y] 


^n-|A| ^ l0g(w - |A|) 
PdA<= PdA'= 


This completes the proof. □ 


Proofs from Section 2.3.2. 

Proof of Theorem 4. We analyze a CTMC Z{t) that 
is stochastically dominated by the total number of infected 
nodes N{t) in the SlSbs epidemic. This being true, we have 
that the expected time to absorption for Z{t) is less than 
that of N{t). We note that this proof structure is similar to 
that of Theorem 4.1 from [8]; however it is much simplified 
via the ergodic embedding technique of [10]. 

Note that for any set of nodes A C U of size |A| < m (for 
some m = n“), we have that the neighborhood of A has size 
at least rim{G)\A\. Now suppose the set of infected nodes I 
satishes j/j < m. Under bounded susceptibility, each node 
in r(/) is infected at a rate which is > 4>(/3). Using this 
observation, we dehne Z{t) to have transition rates given by 

jz —> z -h 1 at rate Pm{G)^{/3)z Vl< 2 <m — 1, 
z ^ z — 1 at rate 2 V1 < 2 < m, 

and all other transitions equal to 0. Note that Z{t) is a finite 
state space CTMC on state space {0,1,... ,m} - standard 


coupling arguments (for example, see [ 12 ]) show that it is 
dominated by N{t). Now using Lemma 1, we have 

m—l 

E[re] > E[inf{t : Z{t) = 0}] = ^ (r,)/(G)$(/l))\ 

k = 0 

So, if ${l3)rim{G) > 1, we have log(E[Te]) = fl(n“). □ 


Now we show that for the Erdos-Renyi graph in the con¬ 
nected regime, the condition for persistent SISba epidemics 
is satisfied with high probability. Note that this proof also 
immediately gives persistence for the SIS epidemic. 

Proof OF Theorem 5. Let G be a realization of a G(n,p) 
graph. Fix a G (0,1),m = n“, and choosing some 7 > 1, we 
dehne b = 7/<l?(/3). We now prove that P [vim{G) < 6 ] < 1/n^ 
for p = i 6 j£gii Then by stochastic dominance, the state¬ 
ment is true for any p > First we have 

m 

^[vl{G)<b]<J2 E P[|r(5)]<H. 

i=l S:|S|=i 


Dehning B{n,p,i) ~ Binomial{n — i, 1 — (1 —pf), we can 
write 


Wm{G) < ^1 < E ( ^]'^[B{n,P,i) < bi] . 


— {n — i)ip{l — ip/2)S^ 


Note that (1—p)* < 1 — ip + Pp'^/2yi. This implies that 
lE[fJ] ^ {n — i)ip (1 — fp/2). 

Dehning 5 = (l - (^„-i)ip{i-ip/ 2 ) ) and using Chernoh bound, 
P [B{n,p, i) < 6i] < P [B{n,p, i) < {1 — S)¥,[B{n,p, i)]] 

<c.Apf -^lB{n,p,i)]5^ ' 

< exp 

Note that (1 — ip/2)S^ —>■ 1 as n —>■ 00 . So for n large 
enough, {n — i)(l — ip/2)S^ > n/2 which implies 

P ['7m(G) < &] < E “ ip/2)S^'j 

m i 

< E ^ exp(-i(np/4)) 

i=l 

m 1 ^ 

= E2Texp(-i(np/4-logn)) < 

• ^ 1 1 


□ 


Next we turn to Theorem 6 . The high-level proof struc¬ 
ture is similar to existing analysis for the SIS epidemics on 
stars [ 8 , 9]. However, to study the effect of bounded suscep¬ 
tibility, we need to carefully characterize various inter-epoch 
events, which require controlling the size of the core and the 
sparsity of the periphery. 

We need two crucial lemmas for the proof. The hrst 
lemma is a crude but general upper bound for SISbs extinction¬ 
time in all graphs, 
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Lemma 4. For any network of size n, and for any arbi¬ 
trary initial set of infected nodes, the expected extinction¬ 
time for SISbs dynamics has an upper bound given by 


E[Te] < 


(1 + 

n^rnax 


Proof. We use the expression for E[re] from Lemma 1. 
For any susceptible node with an infected neighbor, we have 
that the infection rate is bounded by <^max- Further, at any 
time when the number of infected nodes |/| > 1, we can 
stochastically upper bound the infection by assuming that 
all susceptible nodes have an infected neighbor. Thus 


E[re] < 1 +n 


i + 1 


■E 


n^max \ k , 


< 


jl + ^maxT 

n^rnax 


□ 


We use Lemma 4 to control extinction-time in dense sub¬ 
graphs. In case of sparse subgraphs, we can use a result 
from [8], which gives conditions for fast extinction of the 
SIS epidemic (and hence SISf,s epidemic) in terms of p{G). 


Lemma 5. [8, Theorem 3.1] Let p{A) be the spectral ra¬ 
dius of the adjacency matrix A. If P < l/p(A), then the 
probability that the epidemic has not died out by time t, given 
the initial condition X(0), admits the following upper bound. 

F[X{t) / 0] < 

where ||X(0)||i = ^i(O). 


With these two ingredients, we can now prove Theorem 6. 


Proof of Theorem 6. We define cyclical epochs, start¬ 
ing with some (arbitrary) set of nodes in the core 14 infected, 
transitioning to all nodes in 14 becoming susceptible, and 
then finally at least one node in 14 getting reinfected. This 
technique is similar those used in [8, 9] - however instead of a 
star, we now have to deal with core and periphery subgraphs 
and also bounded susceptibility. 

Consider epoch j (for any j G N) starting at time Tqj. 
We define Tij := inf{t > Tqj : '^i^Y^Xi{t) = 0}, i.e., all 
nodes in 14 become uninfected. From Lemma 4, we have 


nTij 


To,j] < 


(1 + ^max)^ 

m^jnax 


( 1 ) 


where m = |f4|. Note that, under bounded susceptibility, 
this only depends on m and not on core/periphery topologies. 

Next, we define Tb.j+i to be the first time after Tij when 
either some node in 14 gets re-infected or all the nodes in Vp 
recover - in case of the former, the next epoch starts, while 
in case of the latter, the infection dies out (i.e., Toj+i = Te). 

It is clear that the random variable Toj+i—Tij is stochas¬ 
tically dominated by 4i, the time taken by all the nodes in Vp 
to completely recover if the periphery was isolated from the 
core at Tij. Now given that the infection profile $ is con¬ 
cave, and $(0) = 0, we have 4?(/fa;) < 4?(/3)a; Va: > 1 and 
so the number of infected nodes under the SlSbs dynamics 
is stochastically dominated by that in an SIS infection pro¬ 
cess with rate ^{P). Now given that ^{P)p{G{Vp)) = 1 — e 
for some e > 0, where p{G{Vp)) is the spectral radius of the 


periphery subgraph, we have from Lemma 5 that P(4i > 
t) < min |ne“ Thus 

^OO 

E[24,j-i-i — Tij] < E[4i] = / P(Ci > t)dt 
Jo 


< 


f 

Jo 




l-HP)p{G{Vp))- 


( 2 ) 


We now derive a lower bound for the probability q that 
the infection dies in the periphery G(f4) before any node 
in the core gets reinfected. Since each node in the core V4 
is infected at a rate at most ^max, the time after Tij for 
the core to get infected can be bounded from below by ( 2 , 
an Exponential(m^rnax) random variable, where m = |f4|. 
Defining to = ^ave 


poo 

IPiCl > C2I = rn^max / P(Ci > 

Jo 


^dt 


< m<f?„ 


-I- m<f>„ 


r 

Vo 

-f 


= 1 — e 
nm^n 


ne 
*0 


-(l-«(/3)p(G(Vp))-|-m«max)t 


dt 


^-(l-'I>(^)p(G(Vp))-|-m'I>max)to 


1 - $(/3)p(G(Fp))-h 
1 - ^{p)p{G{Vp)) 


= 1 - 
= f - 


1 - $(/3)p(G(f4))-h 

e 


— m^rnax 


yp) 

— m^ri 


£ -\- m^max 

The probability q that the infection dies before any node in 
the core gets infected can thus be bounded from below. 

q > P[Ci < C2] > , ( 3 ) 

where G = 1 -I- m^max/e > 1. 

Finally, the extinction-time can now be bounded as 

v-i-i 

Te <st Ui, 


where P is a geometric random variable with mean i and 
Ui are i.i.d. random variables independent of V with mean 
equal to the sum of upper bounds (R.H.S) in Equations f 
and 2. Thus, using Equations 3,1, 2, we have 


E[re] < -E[Ui] 


< Gn 


'^rnax / (1 -f ^maxY 


■ -b 


m^m.ax ' 1 - ^{P)p{G{Vp)) ) ■ 
Simplifying we get that logE[re] = O (mlogn). □ 


4. CONCLUSION 

We have studied the effect of bounded susceptibility on the 
SI and SIS epidemics. In particular, we showed that im¬ 
posing this constraint brings out a sharp dichotomy between 
star-like and clique-like networks. Our results are applicable 
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in several real-life settings which exhibit bounded suscepti¬ 
bility - social networks, human-disease networks, infrastruc¬ 
ture networks, etc. Along the way, we also develop a suite 
of lemmas relating graph structure and epidemic behavior - 
some extend/simplify existing results while others are novel. 
We hope these will prove useful in other settings. 
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APPENDIX 


In this section, we give concentration results which show that 
there is a sharp contrast in spread/extinction times between 
the bounded susceptibility and the traditional models. 

First however, we give a proof of Lemma 1 from Sec¬ 
tion 2.2. This proof is very similar to that of Lemma 8 in 
[10] - we include it here mainly for the sake of completeness. 

Proof of Lemma 1. Following [10], we embed this pro¬ 
cess in another CTMC with same transition rates ex¬ 

cept for an additional transition from 0 to 1 at rate qo,i- It is 
easy to see that {U{t)} is now an ergodic CTMC. Recall that 
we defined T(i,j) := inf{t : U{t) = i]]tf(0) = i}. Similarly, 
dehne f{i,j) := inf{t : U{t) = j][I7(0) = i}. Note that 
r(l,0) and r(l,0) have the same distribution, and there¬ 
fore, E[r(l,0)] = E[T'(1,0)]. Now, {17(t)} is a hnite state, 
irreducible CTMC and therefore, positive recurrent. Let tt 
be its stationary distribution. Solving for ttq gives: 

± = i + y rr^^- 

From standard ergodic Markov Chain theory ]29], we have 
that E[T'(0, 0)] = , and further, from our construction, 

E[f’(0,0)] = l/go,i +E[f'(l,0)]. Thus: 


E[T(1,0)] = ;^ 
< 10,1 



1 

9i.o 


1 +En 


Qi,i + 1 
Qi+l,i 


This completes the proof. □ 


Sha-' Lower Bound 


Theorem 8. For any graph under Shs dynamics, Ta > 
w.p. at least 1- 1=. 

Proof. Since any susceptible node is infected by rate at 
most $max, Tk+i—Tk >at Zk, where Zk has Exponential{{n— 
k)^max) distribution. Using Chernoff bound. 


y^,Zk<t 


< e^ 


'U^[e 

k^l 

n—1 j 

n 


fc -I-1 


Therefore, 


Ta < 


log n 


2$. 


< 




logn 

95 



□ 


SI Shat Lower Bound in Star-like Networks. Before we 
prove that the infection duration does not scale exponen¬ 
tially in star-like networks with high probability, we need 
the following lemma: 

Lemma 6. In a clique of size n, for an arbitrary initial in¬ 
fected node, the second moment of extinction time for SISba 
dynamics has an upper bound given by 


SI: Upper Bound in Clique-like Networks 


E[t,^] < 


2(1 -b max(l, 


Theorem 7. Under the conditions of Theorem 1, the spreading¬ 
time Ta < with probability at least 1 — N 

Proof. Define A(A:) = „ -rr• Then we have: 

'' 7 p 77(G) min(fc,n —fc) 

P [Tfc+i -Tk> A{k)] 

= ^ P[rfc+i-ri,>A(fc)]S'i, = si,]P[S'fc = Sfc] 


Proof. We again consider the evolution of the total num¬ 
ber of infected nodes with transition rates qi,i+i < {n — 
i)^max and qi+i,i = i -I- 1. The second moment for the first 
passage time to 0 starting from state 1 for a general birth 
and death process is derived in [30] and is given by: 


E[r,^] = 


2 

qi.oTo 


E 


1 


Qi—l,i7Li—l 




= E 


exp - 


2p\S{Sk)\logn 


< exp(—2 log n). 


PrjlG) min(fc, n — k) J 
By telescoping, we have Ta = Tfe+i ~ Tfc- Further: 


E; 


qi.OTO qi,i-l 
1 = 1 


1 +En 


<ll,l-ll 
■: 11+1,1 


E^W = E 

fc=l fe=l 


2 logn 

l3rj{G) min(fc, n 


_ < 4(log nf 
k) - prj{G) 


Putting all the above together, we get: 

'n — 1 n — 1 

71-1-1 — Tk > A(fc) 

.fc=l k = l 

n-1 

< Vp[rfc+i-rfc > A(fc)] < -. 
^' n 

k=l 


Ta > 


4(logn)^ 


PviG) 


< 


Using this equation, and qi,i+i < (n — i)$max and qi+i,i = 
i + 1, we have 


E]Te] < 


n^rnaxTO 
2 


< 


E 

i=l 
(1 + ^ 


tvoK 




k fn\ 
max 




71$ 




~)2n 




if ^max > 1 
if ^max < 1 


□ 


This completes the proof. □ 

Therefore, in clique-like networks with ri(G) = co ((logn)^), 
the spread time under the SI model is vanishingly small 
with high probability. Again, as in Proposition 1, we show 
that this is not the case with the Sha epidemic. 


Theorem 9. Under the conditions of Theorem 6, logTe = 

1 — n-.T 

O(m logn) with probability at least 1 — 2n~ —n ’n for 
large enough n. 

Proof. As in the proof of Theorem 6, we consider cycles 
of epidemic evolution. We first bound the total number 
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of cycles and then give a bound for the duration of each 
cycle. As shown in Theorem 6, the total number of cycles is 
stochastically dominated by V, a geometric random variable 
with mean i. Let K = Cn“logn. Therefore, from 
Equation 3, 

P[E > K] = (1 - q)^ < expf-gA") < - (4) 

n 

We next bound the duration of the first epoch in the cycle 
- the duration of the cycle in which at least one node in the 
core is infected. Since this time is stochastically dominated 
by the duration of infection in a clique of size m (size of the 
core), we have 




Using Markov inequality, 
¥[Tij-To,j > 2n 


— A 2m ^ 

4n « 


Using union bound, the total duration of these first epochs 
in K cycles can bounded as follows: 




< K 


2m$2 ji2p. 


< n 2. 


'^max' 
~ ^max 


(5) 


for n large enough. Similarly, we upper bound the duration 
of the second epoch which is stochastically dominated by ((i 
defined in Section 3. 

P[roj+i - Ti,j > ^ logn] < 


Therefore, the total duration of the second epochs in K cy¬ 
cles is bounded by. 


P[^ Toj+i - Tij > ?^logn] < 
i=i 

C log n 1 
nP ^ n 


( 6 ) 


for large enough n. Since Te = Z]J=i To.i+i-Tb.j = Y^]=i{To,j+i- 
Tij) -f {Tij — Toj), using union bound we have, 

P[Te > -f Clogn)] < P[U > K] 

K 

+ ¥[J2Ti,j - To,j > 2An"?*”“*"] 

^ opCC' 

+ - Ti,j > logn]. 

j=i 


Finally, using Equations 4, 5, 6, we have for n large enough, 

P]re > 2A(n~*'"“”’ -|- — logn)] < 2n~^ -|- n 2 ; . 

Since log(A(n~*'"“* -|- ^ logn)) = 0(mlogn), we have our 
result. □ 
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